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GENERATION OF DIFFUSE NON THERMAL ATMOSPHERIC PLASMAS 
Field of the Invention 

This invention relates to the generation of gas plasmas, and more particularly to the 
15 electrically-efficient production of diffuse non-thermal atmospheric gas plasmas. 

Review of the Art Known to the Applicants) 

Nonthermal gas discharges generated at atmospheric pressure find widespread use for 
ozone production (see B. Eliasson and U. Kogelschatz, <fi Nonequilibrium volume 

20 plasma chemical-processing", IEEE Trans. Plasma Science, vol.19, pp. 1063 - 1077, 
1991), pollution control (see B. M. Penetrante, J. N. Barsley, and M. C. Hsaio, "Kinetic 
analysis of non-thermal plasmas used for pollution control", Japan. J. Appl Phys. Vol. 
36, pp. 5007 - 5017, 1997.), and surface modification of Polymer films (see J. 
Friedrich, L. Wigant, W. Unger, A. Lippitz, and H. Wittrich, "Corona, spark and 

25 combined UV and ozone modification of polymer films WeBP23", Surface & Coatings 
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Technology, vol. 98, pp. 879 - 885, 1998.). Typically, such atmospheric plasmas 
consist of many nanosecond-scale microdischarges of streamer-like filaments (see A. 
C. Gentile and M. J. Kushner, "Reaction chemistry and optimization of plasma 
remediation of N x O y from gas streams", /. Appl. Phys. vol.78, pp. 2074 - 2085, 1995.) 

5 and much is known of their fundamental properties and their applications (see also M. 
Baeva, A. Dogan, J. Ehlbeck, A. Pott, and J. Uhlenbusch, "CARS diagnostic and 
modeling of a dielectric barrier discharge", Plasma Chem. Plasma Processing, vol. 19, 
pp. 445 - 466, 1999 and R. Hackam and H. Akiyama, "Air pollution control by 
electrical discharges", IEEE Trans. Dielectrics and Electrical Insulation, vol.7, pp.654 

10 — 683, 2000.). More recently it has been suggested that nonthermal atmospheric gas 
discharges can also be diffuse and luminous with the duration of their discharge current 
pulses in excess of sub-milliseconds (see E. E. Kunhardt, "Generation of large-volume, 
atmospheric-pressure, nonequilibrium plasmas", IEEE Trans. Plasma Sci. 9 vol.28, 
pp.189 - 200, 2000; M. J. Shenton and G. C. Stevens, "Surface modification of 

15 polymer surfaces: atmospheric plasma versus vacuum plasma treatments", /. Physics. 
D: Appl Phys. 9 vol.34, pp. 2761 - 2768, 2001 and F. Massines, A. Rabehi, Ph 
Decomps, R. B. Gadri, P. Segur, and C. Mayoux, "Experimental and theoretical study 
of a glow discharge at atmospheric pressure controlled by dielectric barrier", Journal of 
Applied Physics, vol.83, pp. 2950 - 2957, 1998.). Compared to their counterparts in the 

20 streamer-dominated mode, these diffuse nonthermal atmospheric plasmas have greater 
spatial uniformity, better temporal stability, and much lower gas temperature, typically 
in the range 75°C - 150°C (see A. Schutze A, J. Y. Jeong, S. E. Babayan, J. Park, and 
G. S. Selwyn, 'The atmospheric-pressure plasma jet: A review and comparison to other 
plasma sources", IEEE Trans Plasma Sci, vol. 26, pp.1685 - 1694, 1998.). These 

25 properties make them particularly attractive for a number of key materials processing 
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applications such as etching, deposition, and structural modification of polymeric 
surfaces. These diffuse nonthermal atmospheric plasmas have been generated with 
dielectrically insulated electrodes at audio frequencies (see Massines et al, supra), with 
un-insulated electrodes at radio frequencies (see Schutze et al 9 supra), and at 

5 microwave frequencies (see Shenton and Stevens, supra). Apart from their desirable 
spatial and temporal characteristics, diffuse nonthermal atmospheric plasmas have 
different electrical and chemical properties from that of their streamer-dominated 
counterparts. While many of their fundamental properties remain to be fully 
understood, the focus of research has been largely experimental improvement of 

10 existing applications and empirical exploration of new applications. 

For controlled and optimised performance of material processing, it is always 
desirable to be able to adjust and tailor plasma properties. In the case of diffuse 
nonthermal atmospheric plasmas, this has been pursued mainly through plasma rig 
designs and gas composition, and to a lesser extent through an appropriate choice of 
15 excitation frequency (See J. Park, L Henins, W. Hermann, and G. S. Selwyn, "Gas 
breakdown in an atmospheric pressure radio-frequency capacitive plasma source", J. 
Appl Phys., vol. 89 15 - 19, 2001). 

US Patent number US6228330 teaches the use of atmospheric-pressure plasmas for 
the decontamination and sterilization of sensitive equipment and material, and provides 
20 evidence of their efficacy in inactivation of a number of bacterial spores. The patent 
describes a system for plasma recirculation, but discloses no shaping of the waveform 
of the applied voltage. 
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The European patent application number EP 1 040 839 descibes a multi-stage 
process for sterilization using a gas plasma. The application discloses no shaping of the 
waveform of the applied voltage. 

US Patent number 6,118,218 describes an atmospheric pressure plasma treater 
5 incorporating a porous metallic electrode. The application discloses no shaping of the 
waveform of the applied voltage. 

A number of features of non-thermal gas plasmas are particularly important in 
regard to their applicability. First, the temperature of the gas plasma may be deleterious 
to objects with which the plasma comes into contact, for example during sterilization or 
10 surface decontamination. Second, the power dissipation in the plasma has a direct 
bearing on the energy efficiency of plasma generation, which is important for the 
overall economic efficiency of devices using the plasma, and is particularly important 
in applications remote from a readily-available supply of electrical power. 

Accordingly, it is an object of the present invention to provide a means to generate 
15 atmospheric-pressure plasmas at increased energy efficiency and at low temperature, 
whilst retaining the desirable feature of a high density of active species. 



Summary of the Invention 

In accordance with the inventive concept, these objects are addressed by the use of an 
20 applied voltage that exhibits a waveform (as defined herein) which is truncated (as 
herein defined) and/or which decays asymmetrically from its peak value. 
Advantageously, the applied voltage, V, as a function of time, t, said time t being 
measured from any arbitrary instant, takes the form of a waveform, V(t), of cycle time 
T, wherein in at least one of the half cycles, i.e. . between (t=i T) and (t=i T + T/2) or 
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between (t=i T + T/2) and (1=(i+l)T), where i takes integer values, the waveform is 

characterised by the magnitude of the integral of the voltage with respect to time being 

greater in the first half of said half cycle than in the second half of said half cycle. 

Advantageously also, the applied voltage, V, as a function of time, t, said time t being 
5 measured from any arbitrary instant, takes the form of a waveform, V(t), of cycle time 

T, wherein in at least one of the half cycles, ie. between (t=i T) and (t=i T + T/2) or 

between (t=i T + T/2) and (t=(H-l)T), where i takes integer values, the waveform is 

characterised by a period of substantially constant voltage. 

Also advantageously, the applied voltage may be defined by equation El, below. 
10 Also advantageously, the applied voltage may be defined by equation E2, below. 

Also advantageously, the applied voltage may be defined by equation E3, below. 

Preferably, the applied voltage may be generated by the action of a control system, said 

control system using a measurement of the plasma discharge current as an input signal. 

Included within the scope of the invention is a method of generating a gas plasma 
15 substantially as described herein with reference to and as illustrated by the 

accompanying drawings. 



Brief description of the Drawings 

Figure 1 is a graph showing the reduced ionization coefficient of argon as computed 
20 with a Boltzmann solver (circles) and with the source term technique (4-) together with 
reduced ionization coefficient of nitrogen using the source term technique (solid line). 

Figure 2 is a graph showing typical voltage-current characteristics of a helium-nitrogen 
discharge under sinusoidal excitation with the gas voltage in solid curve, the discharge 
current in thick dashed curve, and the applied voltage in dot curve. 
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Figure 3 is a graph showing a generalised peak-levelled waveform. 

Figure 4 is a graph showing (a) A pulsed excitation voltage (solid curve) at 10kHz 
repetition frequency, made from a sinusoidal voltage (dashed curve) having its peaks 
levelled; and (b) time dependence of the peak-levelled excitation voltage (dot curve), 
5 the gas voltage (solid curve), and the discharge current (thick dashed curve). 

Figure 5 is graph showing normalized plasma power density (circles) and normalized 
electron density (diamonds) as a function of reduced magnitude of the applied voltage 
with 0.5% nitrogen. 

Figure 6 is a graph showing the normalized plasma power density (circles) and 
10 normalized electron density (diamonds) as a function of reduced magnitude of the 
applied voltage with no nitrogen impurities. 

Figure 7 is a graph showing the normalized electron density (circles) and normalized 
metastable density (diamonds) as a function of reduced magnitude of the applied 
voltage with 0.5% nitrogen impurities. 

15 Figure 8 is a graph showing a generalised tail-trimmed waveform. 

Figure 9 is a graph showing (a) a peak-levelled excitation voltage at 10kHz repetition 
frequency (dashed curve) and with its pulse tail trimmed (solid curve); and (b) time 
dependence of the peak-levelled and tail-trimmed excitation voltage (dot curve), the gas 
voltage (solid curve), and the discharge current (thick dashed curve). 

20 Figure 10 is a graph showing (a) a sinusoidal excitation voltage at 10kHz repetition 
frequency (dashed curve) and with its tail shaped with a Gaussian decay (solid curve); 
and (b) time dependence of the excitation voltage (dot curve), the gas voltage (solid 
curve), and the discharge current (thick dashed curve). 
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The Inventive Concept 

The inventive concept considers pulsing plasma excitation voltage as a way to 
control and improve properties of diffuse nonthermal atmospheric plasmas. For 
5 streamer-dominated atmospheric gas discharges, pulsed plasma generation has been 
known to facilitate better energy efficiency and greater control of the glow-to-arc 
transition although this has never been achieved for diffuse, nonthermal atmosphere 
plasmas. If similar improvements could be achieved for such plasmas, it would be 
useful for applications where plasma power consumption is an important issue, for 
10 example aircraft cloaking and industry-scale surface treatment. 

To illustrate the benefits of plasma pulsing in different plasma generation 
configurations, we will consider nonthermal atmospheric gas discharges generated in a 
helium-nitrogen mixture and between two dielectrically insulated parallel-plate 
electrodes. Specifically the inventive concept attempts to answer (1) whether pulsed 

15 plasma generation can lead to saving in electrical power needed to sustain the generated 
plasma; and (2) if so how pulse shape and width may be adjusted to enhance energy 
efficiency. The approach uses on a one-dimensional computer code (See X. T. Deng 
and M. G. Kong, "Parametric conditions for generation of stable atmospheric pressure 
nonthermal plasmas 5 ', in the 2001 Annual Report, Conference on Electrical Insulation 

20 and Dielectric Phenomena 677 (IEEE Catalogue Number: 01CH37225), pp. 677 - 680, 
2001) developed using a hydrodynamic model similar to those employed in most 
numerical studies for diffuse nonthermal atmospheric plasmas (See Massines et al, 
supra and F. Tochikubo, T. Chiba, and T. Watanabe, "Structure of low-frequency 
helium glow discharge at atmospheric pressure between parallel plate dielectric 
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electrodes", Japanese Journal of Applied Physics, vol. 38, pp. 5244 - 5250, 1999). The 
key features of the underlying model of diffuse nonthermal atmospheric plasmas as 
well as its numerical implementation will be described below, for clarity. Voltage- 
current characteristics of nonthermal atmospheric plasma under sinusoidal excitation 
5 will be used as an example to demonstrate the utility of this approach. 



A physical model of nonthermal atmospheric plasmas 

Diffuse nonthermal atmospheric plasmas may be induced and sustained between two 
parallel-plate electrodes, each one optionally coated with a dielectric layer and 

10 connected externally to a sinusoidal voltage source, typically at voltages in excess of 
lkV, and at audio frequencies. The background gas is, for example, atmospheric 
helium mixed with a small fraction of nitrogen (up to 1%) at room temperature of 
293K, and the dynamics of the generated nonthermal atmospheric plasma is described 
by the Boltzmann equation, well known in the art It has been established that the 

15 hydrodynamic assumptions can be applied to diffuse nonthermal atmospheric plasmas 
(See Massines et al and Tochikubo et al, supra). This reduces the Boltzmann equations 
to the continuity and momentum transfer equations for electrons and ions, thus 
facilitating a macroscopic description of plasma dynamics without its microscopic 
details. On the other hand, the dynamics of electrons and ions are determined by the 

20 electric field in the space between the two parallel electrodes, which has two 
components, one induced by the externally applied excitation voltage and the other by 
space charges. The electric field can be calculated by solving Poisson's equation. In the 
one-dimensional limit, the Boltzmann equations and Poisson's equations are closely 
coupled together as follows 
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(la) 



an, _ c/ A djn^tWAnt)] d 2 [n,{r 9 t)D t {r 9 t)] 
— -S, (r, t) - + Qr 2 



(lb) 




(lc) 



where n+ and n. are the ion and electron densities respectively, and m is the density of 
5 the z ,th neutral species considered. S± and Si are the source terms for charged particles 
and neutral species respectively. D and W are diffusion coefficient and drift velocity 
with subscripts -, +, and i denoting respectively electrons, ions, and i th neutral species 
considered in the physical model. E is the electric field, \e\ is the charge of the electron, 
s 0 is the dielectric permittivity, and r is the variable representing the spatial position 
10 normal to the electrode. With the hydrodynamic assumptions, reaction rate coefficients, 
ionization coefficients, drift velocity, and diffusion coefficients can be approximated as 
a function of the electric field in the ionised gas between the two electrodes. Thus this 
physical model is similar to those used in Massines et al {supra) and Tochikubo et al, 
(supra). The values of these coefficients are obtained from relevant experiments and 
15 several rate compilation studies which will be discussed in greater detail below. These 
allow eqs.(la) and (lb) to be solved to yield densities of electrons, ionic species, and 
metastables, which then allows through eq.(lc) calculation of electric field. Boundary 
conditions between the ionized gas and the dielectrically coated electrodes need to be 
carefully treated To this end, a circuit equation is included to relate the electric field in 
20 the gas to the source voltage (the output voltage of the power supply) via a source 
resistor, R S9 and two serial capacitors each representing one dielectric coating layer. The 
numerical algorithm used to solve the above equations is essentially based on the 
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Patankar scheme, and the discretization employs the upwind scheme (see S. Patankar, 
"Numerical heat transfer and fluid flow". Hemisphere publishing Co, 1980). 

The model considers reactions involving eight different species, namely (1) two 
ground-state neutral species, He(l ! S) and N2(4°S); (2) two helium metastables, He(2 3 S) 

5 and He(2 1 S); (3) electrons; (4) ground-state atomic helium ions, He + ; (5) ground-state 
molecular helium ions, He2 + (2 3 SS U + ); and (6) ground state molecular nitrogen ions, 
N2 + . In total, the model considers 17 reactions, as detailed in Table 1 below, together 
with their reaction rates and reference sources. Products of these 17 reactions include 
two additional species, namely excited molecular helium (He*) and atomic nitrogen 

10 (N). Any subsequent reaction of these two additional species with the original eight 
species is not considered in the model. It is worth noting that the model is capable of 
evaluating production of molecular nitrogen ions that are not considered in reported 
numerical studies of dififiise nonthermal atmospheric discharges (e.g. Massines et al 9 
supra and Tochikubo et al 9 supra). 

15 



Table 1: Reactions considered for a He-N2 discharge and their rate coefficients 



No 


Reaction 


Reaction rate 


Reference 


Direct ionization 


1 


He + e->He + + e + e 




1, 2,3 


2 


N 2 + e N 2 + + e + e 


a 2 


1,4 


Excitation 


3 


He + e -> He (2*S) + e 


A 


5 


4 


He + e -» He (2*S) + e 




5 


5 


He (23S) + 2He (1*S) -> He 2 * + He (VS) 


1.9 X 10- 34 cm 6 /s 


6 


De-excitation 


6 


He (23S) + e -» He (l l S) + e 


4.2 XlO^ cm 3 /s | 7 


Penning ionisation 


7 


| He (23S) + N 2 He (VS) +N 2 + + e 


8 XlO 11 cmVs 


5,8 


Stepwise ionisation 
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No 


Reaction 


Reaction rate 


Reference 


8 


He (23S) + He (2*S) -> He (l^S) + He + + e 


2.9 XlO" 9 cm 3 /s 


9 


Charge transfer 


9 


He + + 2He (l'S) -» He 2 + + He (l'S) 


6.3 XlO" 32 cmVs 


13,5 


Recombination 


10 


He + + e He (1'S) 


2 XlO" 12 cmVs 


10 


11 


He 2 + + e -» He 2 * 


5 XlO" 16 cm 3 /s 


11 


12 


He+ + e + He He (1«S) + He (2»S) 


1 XlO" 27 cm 6 /s 


10 


13 


He2++2e->He(l*S) + e 


7.1 XlO" 20 cm 6 /s 


13,5 


14 


He 2 + + e -» 2He (23S) + He 


5 XlO' 9 cmVs 


6 


15 


He 2 + + e +He -» He 2 * + He 


5 XlO** 7 cm 6 /s 


11 


16 


N 2 + + e -> N + N 


4.8 XlO" 8 cmVs 


8 


17 


N 2 + + 2e -» N 2 + e 


1.4 XlO" 26 cm 6 /s 


12 



10 



15 
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(3) J. Dutton, "A survey of electron swarm data", J. Phys Chem. Ref. Data, vol.4, 
pp. 577-856,1975. 

(4) I. Peres, N. Ouadoudi, L. C. Pitchford, and J. P. Boeuf, "Analytical 
formulation of ionization source term for discharge models in argon, helium, 
nitrogen, and silane", J. Appl. Phys., vol.72, pp. 4533 - 4536, 1992. 

(5) R. Ben Gadri, ''Numerical simulation of an atmospheric pressure and dielectric 
barrier controlled glow discharge", PhD thesis, University Paul Sabatier of 
Toulouse, France, 1997. 

(6) J. M. Pouvesle, A. Bouchoule, and J. Stevefelt, "Modeling of the charge 
transfer afterglow excited by intense electrical discharges in high pressure 
helium nitrogen mixtures", J. Chem. Phys., vol.77, pp. 817 - 820, 1982. 
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(7) R. Deloche, P. Monchicourt, M. Cheret, and F. Lambert, "High-pressure 
helium afterglow at room temperature", Phys. Rev. A, vol. 13, pp. 1140 - 
1176, 1976. 

(8) A. Cenian, A. Chernukho, and V. Borodin, "Modeling of plasma-chemical 
5 reactions in gas mixture of C02 laser", Contribution to Plasma Physics, 

vol.35, pp. 273 -296, 1995. 

(9) L. L. Alves, G. Gousset, and C. M, Ferreira, "A collisional-radiative model for 
microwave discharges in helium at low and intermediate pressures", J. Phys. 
D: Appl. Phys., vol. 25, pp. 1713 -1732, 1992 

(10) T. Quinteros, H. Gao, D. R. DeWitt, R. Schuch, S. Pajek, S. Asp, and D. 
Belkic, "Recombination of D+ and He+ ions with low energy free electrons", 
Phys. Rev. A, vol.51, pp. 1340-1346, 1995. 

(1 1) P. C. Hill and P. R. Herman, "Reaction processes in a He2 + ( 2 II U a 2 S g + ) flash 
lamp", Phys. Rev. A., vol.47, pp. 4837-4844, 1993. 

(12) I. A. Kossyl, A. Yu Kostinsky, A. A. Matveyev, and V. P. Silakov, "Kinetic 
scheme of the non-equilibrium discharge in nitrogen-oxygen mixtures", 
Plasma Source Science and Tech, vol.1, pp. 207 - 220, 1992. 

(13) Tochikubo etal 9 {supra). 

1 Ionization coefficients in helium can be evaluated using Ward's formula (A. L. Ward, 
supra) or experimental data compiled by Dutton (supra). Both agree well with those used 
in Massines et al (supra) , particularly at large reduced electric field (>50Td). Our 
model uses Ward's formula. On the other hand, it is known from spectroscopic 
measurements that discharge dynamics in atmospheric helium is very strongly affected 
• by the presence of impuries, particularly nitrogen. For diffuse nonthermal atmospheric 
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helium discharges, the presence of nitrogen impuries has been either ignored or 

implemented with its ionization coefficient approximated by that of argon. Since argon 
^ has a greater ionization coefficient that nitrogen, the approximation of using argon 

ionization coefficient overestimates the ionization of nitrogen impurities. To this end, 
5 we note that it is possible to formulate analytically ionization source terms for 

discharges in argon, helium, and nitrogen (See Peres, supra). Thus under the same 

conditions 

a A _ ionization source term for species A 
a B ionization source term for species B 

where a is Townsend's first ionization coefficient with its subscript A and B indicating 
10 respectively species A and B. If species A is argon and species B is helium, the above 
formula can be used to predict the ionization coefficient of argon from that of helium as 
calculated from Ward's formula. As shown in Figure 1, the argon ionization coefficient 
calculated with this technique, indicated by crosses 1 leads to an excellent agreement 
with that computed using a Boltzmann solver (see Massines et al, supra), indicated by 
15 circles 2. Therefore equation 2 offers a simple yet reliable way to estimate the 
ionization coefficient of one gas from the known ionization coefficient of a reference 
gas. Similarly if we assume species A is nitrogen and B is helium, the ionization 
coefficient of nitrogen is obtained and this is shown in Figure 1, indicated by the solid 
line 3, as a function of the reduced electric field. The ionization coefficient thus 
20 calculated is used in our model for nitrogen (reaction #2 in Table 1). For excitation 
coefficients, we use that employed by Ben Gadri (st(prd). Similarly our model employs 
the same data used by Ben Gadri {supra) for diffusion coefficients and drift velocity. 
By way of example, and demonstration of the utility of the model presented above to 
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accurately and correctly predict the evolution and properties of atmospheric diffuse 
nonthermal gas plasmas, a simulation is now presented, as follows. 

In this example, we consider a system with the two electrodes having a common radius 

5 of 2cm and separated with a gap fixed at 0.5cm. The total capacitance of the dielectric 
coatings on the two parallel-plate electrodes is 70pF, and the source resistoT in the 
external circuit is 50£l. In this reference case, we consider a sinusoidal source voltage at 
10kHz with a peak voltage of 1.5kV. Before the power supply is switched onto the 
plasma rig, the two electrodes are charged at 950V and so the initial gas voltage is 

10 950V. The hackground gas is 99.5% helium with 0.5% nitrogen unless otherwise 
stated, and its breakdown voltage is assumed to be l.lkV, similar to the choice in 
Massines et al (supra). The secondary electron emission coefficient is dependent on the 
surface condition of electrodes and as such it is not possible to choose a reliable 
coefficient that is applicable to most cases. Numerically different secondary emission 

15 coefficients in the 0.01 — 0.2 range affect the peak value of the discharge current For 
numerical studies considered here, we choose 0.2 for atomic helium ions, 0.1 for 
molecular helium ions, and 0.01 for nitrogen ions (see A. von Engel, "Ionized Gases'% 
Chapter 3, reprinted by American Institute of Physics, 1994). Our numerical model 
assumes that diffuse nonthermal atmospheric plasma is established if its voltage and 

20 current are continuous and repetitive over at least 10 cycles of the applied voltage 
signal. Noting that discharge currents of both thermal plasmas and streamer-dominated 
nonthermal atmospheric plasmas tend to exhibit microscopic current pulses in tens 
nanosecond or less, pulse duration of the discharge current is an additional indicator for 
diffuse atmospheric plasmas. Figure 2 shows a plot of the applied voltage 4, the gas 

25 voltage 5, and the discharge current 6 as a function of time. Though not shown in 
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Figure 2, numerical results confirm that the discharge current 6 is in fact identical 
through many tens of cycles of the applied voltage 4 and so the discharge plasma is 
stable and repetitive temporally over a long period of time. The voltage-current 
characteristics are very similar to that obtained in a comparable experiment and its 
5 numerical simulation (Massines et al, supra), with a typical pattern of one discharge 
every half cycle of the applied voltage. Both gas voltage 5 and discharge current 6 
exhibit very similar waveforms as those observed by Massines et al (supra) and 
Tochikubo et al (supra). The peak discharge current 6 in Figure 2 is between 52 - 58 
mA, lower than 90mA calculated by Massines et al. As shown in Figure 1, this may be 
10 partly due to the use of the larger ionization coefficient of argon used by Massines et al 
to approximate the nitrogen ionization coefficient. It is interesting to note that the peak 
current density of 7.2mA/cm 2 (=-90mA/7r2 2 ) is the highest reported for diffuse 
atmospheric helium discharges generated between dielectrically coated electrodes and 
at kilohertz frequencies. In a series of very similar experiments performed at Tennessee 
15 University, the peak current density is found at most ImA/cm 2 (J. R. Roth, ''Industrial 
plasma engineering", Chapter 12, vol.1, IoP publishing, Bristol, 1995.). A similar peak 
current density of about 0.4mA/cm 2 (=11.25mA/7r3 2 ) was measured with 3kHz plasma 
excitation by the Okazaki group at Sophia University in Japan (T. Yokoyama, M. 
Kogoma, T. Moriwaki, and S. Okazaki, "The mechanism of the sterilisation of glow 
20 plasmas at atmospheric pressure", J. Phys. D: Appl Phys., vol. 23, pp. 1125 - 1128, 
1990), whereas another Japanese group measured 1.8mA/cm 2 with a 100kHz power 
supply. More recently a helium discharge experiment performed at Minnesota 
University recorded a peak current density measured at 0.5mA/cm 2 with 15kHz driving 
voltage (L. Mangolini, K. Orlov, U. Kortshagen, H. Heberlein, and U. Kogelschatz, 
25 "Radial structure of low-frequency atmospheric-pressure glow discharge in helium", 
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Appl Phys. Lett, vol. 80, pp. 1722 - 1724, 2002). Our computed current density of 
4.4mA/cm 2 (=55mA/^2 2 ) falls between these experimental measurements, indicating a 
good current prediction capability of our model. 

Electric power consumption in the plasma is found to be 298mW/cm 3 from our 

5 model. This is almost identical to 300mW/cm 3 measured by Massines et al (supra) and 
similar to 277mW/cm 3 measured by Chen et al (Z. Chen, J. E. Morrison, R. Ben Gadri, 
and J. R. Roth, "A low-frequency impedance matching circuit for a one atmospheric 
uniform glow discharge plasma reactor", paper 6P63, presented at the 25 th IEEE 
International Conference on Plasma Science, Raleigh, USA, June 1998). We have also 

10 calculated electron density and ion density across the space between the two electrodes. 
The time and spatial variation of electron density is similar to that computed in 
Massines et al (supra). A rapid rise of electron density occurs near 250//S and is a result 
of the gas ionisation (or discharge) in that half cycle (250 - 300/^s), and the relatively 
uncharged electron density peak during the remaining period of the half cycle suggests 

15 a phase of inactive electron production. The peak electron density predicted with our 
model is 0.95xl0 10 cm' 3 , not very dissimilar to 2xl0 10 cm" 3 by Tochikubo et al (supra) 
but much lower than 30xl0 10 cm" 3 calculated by Massines et al (supra). Again the large 
electron density reported by Massines et al (supra) may be a result of their 
overestimated ionization coefficient for nitrogen (see Figure 1). For diffuse atmospheric 

20 helium plasmas, there is very little data of direct measurement of electron density for 
comparison with our calculated electron density. However indirectly the ion-trapping 
mechanism proposed by Roth (supra) may be used to estimate the average electron 
density and this can be captured by the following equation: 

2mv c P„ 

e e 2 E 2 0 K) 
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where n e is the average electron density, m electron mass, e electron charge, v c electron 
collision frequency in helium, P av the average electric power density consumed in the 
plasma, and E 0 the peak electric field in the plasma. Assuming P av =300mW/cm 3 , 
£ 0 =3kV/cm, and v c =1.8xl0 12 Hz (see Roth, supra), eq.(3) yields 4.3xl0 8 cm~ 3 . Given 
5 that the peak electron density is approximately 10-50 times greater than the calculated 
average electron density, the above calculation suggests that the peak electron density is 
between 4xl0 9 cm' 3 and 2xl0 10 cm~ 3 . Thus our calculated peak electron density appears 
reasonable. This is also confirmed by electron density estimate using the current density 
formula of J = en&E (where n is the electron mobility) as employed for measuring 

10 electron density in micro-hollow cathode discharges (see R. H. Stark and K. H. 
Schoenbach, "Direct current high-pressure glow discharges", Journal of Applied 
Physics, vol. 85, pp. 2075 - 2080, 1999). In general our simulated voltage-current 
characteristics, the peak discharge current density, the plasma power density, and the 
electron density yield a favourable comparison with available experimental and 

15 numerical data of comparable diffuse atmospheric helium discharges. While there is 
much scope to include other features of diffuse atmospheric plasmas, for example 
multi-dimensional effects, our plasma model described above is capable of capturing 
the main plasma features accurately and as such it will be used to explore the benefits 
of plasma pulsing. 

20 

The underlying inventive concept, realising a benefit of energy saving is best 
understood from the voltage-current characteristic in Figure 2 (Massines et al, sttprd). It 
is seen that the discharge current 6 is of very low magnitude after its peak between 
207//S - 245//S in the half-cycle from 200//s to 250/zs, when the applied voltage 4 is 
25 very la/ge. This suggests that between the instant of this discharge current pulse and 
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the polarity reversal point (that occurs at around 250/is) of the applied voltage 4 in the 
same half cycle the applied voltage 4 does not necessarily contribute to plasma 
generation but may adversely heat up electrons. Through energy transfer from electrons 
to heavy particles, this undesirable electron heating would increase gas temperature. 
5 Therefore by reducing the applied voltage 4 over this "inactive" period of electron 
production, the gas temperature may be lowered and equally importantly the input 
power can be reduced. 



Description of the Preferred Embodiments 
10 We will now describe preferred embodiments of the invention by means of four 
examples. The first uses a waveform essentially comprising a peak-levelled sinusoid. 
The second employs a 'tail-trimmed' sinusoid, and the third employs a 'tail-shaped* 
sinusoid. The fourth example describes how a control system can be used to shape the 
waveform of the applied voltage in response to the resultant discharge current 

15 

Example 1 -Peak-Levelled Sinusoid 

We first consider a simple pulsed excitation voltage waveform, shown in Figure 
4a, which is essentially derived from a sinusoidal wave with its peak (positive and 
negative) levelled to a flat top. Voltage waveform of Figure 4a can be easily obtained 
20 electronically, and so it represents a realistic option. By reference to Figure 3, the 
general equation describing the peak-levelled pulse waveform, V(t), is as follows: 
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< remain f t, 2k) < 




U 'J ■ 



< rem(2n f t 9 2ri) < 



V m sm(27t ft) 



elsewhere 



where : 




and: 

rem(2 n f t 9 In) gives the remainder after the division (2 n f t) I In . 

5 

With reference to Figure 4, the magnitude of the original sinusoidal voltage, V S9 is fixed 
at 1.5kV, whereas that of the peak-levelled voltage 7, V p9 varies from 0.4kV to 1.5kV. 
The repetition frequency of the excitation voltage remains at 10kHz. 

Figure 4b shows the voltage and current characteristics of a diffuse nonthermal 

10 atmospheric gas discharge generated with V p = 1.2kV. The discharge current 8 remains 
repetitive and stable through many tens of cycles of the applied voltage 7. Therefore 
under the conditions considered here, pulsing the plasma-generating voltage 7 does not 
appear to significantly affect the establishment of diffuse nonthermal atmospheric 
plasma. More specifically the waveform of the discharge current 8 remains relatively 

15 unchanged from that under the sinusoidal excitation, and the peak current is again 
around 55 mA very similar to the sinusoidal case of Figure 2. 

The basic voltage-current characteristics remain relatively unchanged until V v 
decreases below 750V when the generated plasma becomes unsustainable and 
eventually extinguishes. This is illustrated in Figure 5 where the normalized electron 

20 density shown by diamonds 9 and the normalized plasma power density shown by 
circles 10 are plotted against the normalized magnitude of the applied voltage, Vp/V s . 
The normalizations of electron density and plasma power density in this figure, and in 
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Figures 6 and 7 are carried out with respect to the electron density and plasma power 
density obtained by sinusoidal excitation, as illustrated in Figure 2). Lines 11 and 12 
are drawn through the electron density data 9 and the plasma power density data 10 
respectively, to guide the eye. It is evident that when Fp/F s > 0.5 plasma pulsing does 
5 not significantly affect electron density while the plasma power is reduced by up to 
40%. Further numerical simulations suggest that ion densities are also largely 
unaffected as long as Vp/V s > 0.5. Given that electron and ion densities remain 
approximately the same and that the discharge current and the gas voltage undergo 
relatively small changes, the basic plasma characteristics should remain relatively 
10 unchanged. In other words voltage pulsing in Figure 4 is unlikely to affect adversely 
plasma-activated applications that rely on the production of electrons or/and ions, yet is 
capable of considerable power saving. 

It is found that the amount of power saving is dependent on that of nitrogen 
impuries. In the case of very little nitrogen impurities (<0.01%), numerical studies 
15 suggest a power saving of more than 50% as indicated in Figure 6. Again such 
significant power saving is achieved without affecting the production of electrons and 
ions. Also found is that the power saving is not at the expense of the production of 
excited neutral species either. With 0.5% nitrogen, Figure 7 shows the normalized 
density of He(2 3 S) 13 as a function of the V^V % ratio. It is evident that the trend of the 
20 metastable density 13 effectively tracks that of the electron density 14. So if the plasma 
power is reduced with the electron density 14 kept approximately the same as that with 
sinusoidal excitation (without pulsing), the metastable density 13 will also remain at 
approximately the same level as its value with sinusoidal excitation. 

The density tracking between electrons and metastables suggests a key role 
25 played by electrons in influencing production of excited neutral species in diffuse 
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atmospheric gas discharges. It is known that energy required to excite neutral species 
and create chemically reactive species is usually provided through the kinetic energy of 
electrons. There are also reactive species whose densities are directly proportional to 
electron density. For example, densities of most oxygen species (e.g. atomic oxygen, 
5 singlet-sigma metastable oxygen, singlet-delta metastable oxygen) in a diffuse 
nonthermal atmospheric discharge in an oxygen-helium mixture are found to 
proportional to P£? where m = 0.8 - 1.9 and JVis the input RF power to the generated 
plasma (see J. Y. Jeong, J. Park, I. Henins, S. E. Babayan, V. J. Tu, G. S. Selwyn, G. 
Ding, and R. F. Hicks, "Reaction chemistry in the afterglow of an oxygen-helium 

10 atmospheric plasma", J. Phys. Chem., vol. 104, pp. 8027 - 8032, 2000). According to 
eq.(3), these oxygen densities are also proportional to electron density. Therefore the 
very similar electron densities in the sinusoidal case of Figure 2 and the pulsed case of 
Figure 4 may be a good basis to consider that the densities of most reactive species in 
the two cases would be similar and so would their chemical reactivity. Thus it is 

15 possible to produce at least some reactive species electrically efficiently with voltage 
pulsing, and as such this will also be useful for applications that reply on reactive 
species. It is important to note however that in general there is no direct correlation 
between electron density and densities of all reactive species. There are reactive species 
whose densities may be inversely proportional to electron density, for example ozone in 

20 the oxygen-helium plasma mentioned Jeong et al {supra). So for individual diffuse 
nonthermal atmospheric plasmas, the suggested maintenance of production of reactive 
species deduced from that of electron production needs to be assessed carefully. Our 
observation is that it is possible for diffuse nonthermal atmospheric plasmas to maintain 
densities of some reactive species with reduced electric power. 
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Example 2-Tail-Trimmed Waveform 

An alternative embodiment to achieve the benefits of the invention may be 
referred to as a tail-trimmed pulse waveform. With reference to Figure 8, such a 
waveform may be described by the following equations: 



and where ceil(2 f t) rounds (2 / 1) to the nearest integer towards infinity. 

With reference to Figure 9a, the "tail-off' phase of the applied voltage in 
(between 25//s and 50//s) may be trimmed to further enhance energy efficiency for 
plasma generation. The voltage-current characteristics of the induced atmospheric 
10 plasma from this waveform are shown in Figure 9b, where 15 is the applied voltage, 1 6 
is the gas voltage, and 17 is the discharge current It is evident that moderate voltage 
trimming does not significantly affect the generation and characteristics of induced 
nonthermal atmospheric plasmas, though the discharge current 17 has different peak 
values in different half cycles. Numerical studies suggest that further power saving is 




where : 



5 




0 <> x 2 £ 2x x 




WO 2004/016052 PCT/GB2003/003485 

23 

possible but very modest, typically a few percent and at best 10% before electron 
density starts to decrease. 

We have so far established that generation of diffuse nonthermal atmospheric 
plasmas can be made more efficient by altering the waveform of the applied voltage 
5 without affecting the basic plasma characteristics. 

Example 3-Tail-Sliaperi Waveform 

In this example, we consider how pulse width may affect plasma generation. In 

this example we consider one particular type of pulsed voltage signal constructed from 

10 a sinusoidal signal for the voltage rise phase and a Gaussian decay for the voltage tail 

phase. Mathematically in each cycle they may be expressed as follows 

Xsinor if 2Nn < t < t 0 + 2Nn 

r p exp[-(r-< 0 -2j\fo) 2 /T 2 ] if t 0 +2NttZt <2Na + n „ 

V( - t)z= \V 0 sm(Dt if 2Nn + n<.t<t 0 +2Nx + x 

-F p exp[-(*-/ 0 -2iV>r-7z-) 2 /r 2 ] if t 0 +2Nx + x <,t < 2Nn+2n 

where V 0 is the peak voltage of the sinusoidal signal, V v is the peak voltage of the 
Gaussian decay signal, t is the pulse width of the Gaussian signal, a> = 2n/T is the 
15 angular frequency of the applied voltage with T being its period, and t 0 the instant at 
which the sinusoidal and the Gaussian signals joint. T 0 and t 0 are defined as follows: 

t 0 =T 0 +kT/2 (5) 
k = mod{t,T/2) 

To compare to the sinusoidal case of Figure 2, we set Vo = 1.5kV, V p = 1.4kV, and 
r=T/8. The waveform of the applied voltage 18 is shown in Figure 10a, and the voltage 
20 and current signals of the generated atmospheric plasma is shown in Figure 10b, where 
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19 is the applied voltage, 20 is the gas voltage, and 21 is the discharge current It is 
evident that the pattern of one discharge every half cycle remains and this is repetitive 
for many tens of cycles. On the other hand, the gas voltage is markedly different from 
that in the sinusoidal case of Figure 2. Particularly it has a large hump preceding the 
5 main peak that causes the discharge event in each half cycle. Also after the discharge 
event where the discharge current peaks, the gas discharge does not reduce to zero as 
rapidly as in the sinusoidal case. 

The discharge current 21 is between 30 — 45mA, markedly lower than 55mA 
calculated for the sinusoidal case. Interestingly the pulse width of the discharge current 

10 21 is much larger than that in the sinusoidal case (Figure 2), particularly for the positive 
half cycles. For the case in Figure 10, electric power density consumed in the plasma is 
about 0.32W/cm 3 , about 7% above that in the sinusoidal case. Further numerical 
examples studied suggest that under other pulsing conditions diffuse atmospheric 
helium plasmas excited with the pulsed voltage of Figure 10a consume approximately 

15 the same amount of electric power as those with sinusoidal excitation. Therefore 
shortening the pulse width may not be as effective in achieving energy saving on an 
absolute basis, but, as discussed below, have advantages in energy-saving for a given 
density of electrons and metastables achieved. From the voltage-current characteristic, 
this is likely due to the wider pulse widths for both the discharge current and the gas 

20 voltage. On the other hand, electron density achievable can be much higher. For the 
case shown in Figure 10b, the peak electron density is found to be 1.74xl0 l0 cm" 3 , 83% 
above 0.95xl0 10 cm 3 for the sinusoidal case. The spatial variation of electron density 
over one half cycle shows a similar distribution to that that produced by the sinusoidal 
waveform presented above yet with clear difference near the voltage polarity reversal 

25 point. There is a clear evidence of pulsed production of electrons. Further numerical 
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calculations suggest that by changing the pulse width of the applied voltage in Figure 
10a, electron and metastable densities can be made much greater than that under 
sinusoidal excitation. Though inappropriate wave-shaping of the applied voltage can 
also reduce electron and metastable densities, the above discussions suggest that 
5 through control of pulse width it is possible to enhance production of electrons and 
metastables for a given plasma power consumption. To see this more clearly, we 
calculate the ratio of electron density to plasma power density (number of electrons that 
can be produced at a given input power). This is 5.4x1 0 10 /W for this example, and 
3.2xl0 l0 /W for the sinusoidal waveform presented above, an improvement of 68%. 
10 From this comparison, it is still electrically efficient to pulse the applied voltage for the 
generation of diffuse nonthermal atmospheric plasmas. 



Example 4 - Controlled Waveform 

It is clear that an optimal saving of energy in the production of gas plasmas may be 
15 obtained by the use of a control system, such as a feedback control system, wherein the 
waveform of the applied voltage is driven by the control system, which itself uses the 
discharge current as an input variable. In this way, the voltage signal may be attenuated 
following the plasma discharge, so as to make most efficient use of the input energy. 
The use of such a control system can then also take account of any temporal changes in 
20 the system configuration (such as changes in electrode or dielectric coating 
characteristics, or changes in the gas composition), leading to waveforms (as defined 
herein) that may vary in their repetition rate, and which may vary in peak shape and 
amplitude from cycle to cycle. 



25 
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Parametric Ranges of Key System Variables 

Electrode size: Typically between 50cm 2 and 100cm 2 of either circular or rectangular 
shapes. Scaling up is possible, with the electrode area divided into different sections 
connected in parallel to the supply power source. To control plasma stability with large 
5 electrodes, it may be useful to use individually valued series resistors or impedance 
networks in each of these parallel circuit branches. On the other hand, plasma stability 
is easier to achieve with smaller electrodes than 50cm 2 . 

Electrode separation: Typically around 1cm. Smaller electrode separation makes it 
easier to control plasma stability. Although larger electrode separation tends to un- 
10 stabilise the generated nonthermal plasma, it is quite possible to increase the electrode 
separation up to tens of centimetres, especially when the electrode area is large and the 
electric power pulsed. 

Capacitance of dielectric coatings: For the case where at least one of the electrodes has 
a dielectric coating, their capacitance depends on the size of the electrodes and hence 

15 that of the dielectric. Typically for a surface area of 10cm 2 the capacitance of one 
dielectric coating is between lOpF and lOOpF. If the surface area of the dielectric 
coatings is increased, the capacitance increases proportionally. Also it is equally 
acceptable to coat either one electrode or both electrodes, the latter of which will have 
the coating capacitance halved. 

20 Peak voltage: This depends on the repetition frequency, the electrode separation, and 
gas composition. At a repetition frequency of 10kHz and an electrode separation of 
0.5cm, the peak voltage needs to be greater than l.OkV but usually less than 2kV in 
helium or helium mixed with a small fraction of nitrogen or/and oxygen (typically less 
than 1%). This needs to be increased when the electrode separation is increased. At a 

25 greater repetition frequency on the other hand, the peak voltage required becomes 
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smaller. For example between 10 and 20MHz, the peak voltage can be reduced to as 
low as 500V in helium with 0.5cm electrode separation. 

Gas composition: Our numerical work was performed for helium-nitrogen mixture. For 
applications, helium-oxygen is more useful and argon-oxygen is even better as argon is 
5 cheaper than helium. The most practically convenient gas is air (though very difficult to 
simulate). Alternatively noble gases mixed with a small fraction of a reactive gas (such 
as, but not limited to, oxygen can be used. 

to in the Gaussian decay Junction: to is given in equation 4 and 5 in Example 3 of the 
preferred embodiments. Typically it is between 5% and 15% of the repetition period of 
10 the applied electric voltage. 



Summary 

Diffuse nonthermal gas discharges generated at atmospheric pressure have found 
increasing applications in many key materials processing areas such as etching, 

15 deposition, and structural modification of polymeric surfaces. To facilitate tailored 
and improved applications of these novel gas plasmas, we consider their pulsed 
generation based on one-dimensional numerical simulation of helium-nitrogen 
discharges. We consider four waveforms of the plasma-generating voltage, namely 
(1) sinusoidal; (2) peak-levelled sinusoidal; (3) peak-levelled and tail-trimmed 

20 sinusoidal; and (4) pulsed with a Gaussian-shaped tail, all at the same repetition 
frequency of 10kHz. For each case, voltage and current characteristics are calculated 
and then used to assess whether the generated plasma is diffuse and nonthermal. 
Densities of electrons, ions, and metastables are calculated, together with the 
dissipated electric power in the plasma bulk. It is found that plasma pulsing can 

25 significantly reduce the electric power needed to sustain diffuse nonthermal 
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atmospheric plasmas. Specifically by choosing appropriate pulse shape, the plasma- 
sustaining power can be reduced by more than 50% without reducing densities of 
electrons, ions, and metastables. On the other hand, electron density can be 
enhanced by 68% with the same input electric power if the pulse width is suitably 
5 narrowed. 

Scope of the Invention 

The invention is defined in the Claims, which now follow, and in which the term 
"waveform" is understood to include periodic signals that have a zero value over 
10 more than an instantaneous part of the period, and which are more commonly 
referred to as "pulse trains". 

In the Claims, the term "truncated" is taken to include limited in amplitude to a 
maximum and/or minimum voltage, such limitation typically leading to a waveform 
with a substantially flat profile at its extreme value or values, and/or to include 
15 narrowed in pulse width. 

In the Claims, Equation El is taken to be: 




F M sin(2*r/0 



elsewhere 



where : 




and: 



20 



rem(2 7t f t, 2n) gives the remainder after the division (2 n f )ll7t . 
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and further, Equation E2 is taken to be: 



K m sin(|-*,) 

V m sin(| - x, )sin(2 it ft + 0) 

V m sin(?-- x,)sin(2 n ft + &) 
V m sm(2n ft) 



^-x^<rem(2n ft, 2^)£^+x, -* 2 



elsewhere 



where : 



1^2 # /rj 

0 = cetf(2/»)*(l-yJ 

nd where cez/(2 / 0 rounds (2 / r) to the nearest integer towards infinity. 
5 and Equation E3 is taken to be: 



V(t) = 



if 2Nn£t<t 0 +2Nir 

if t 0 +2Nn£t<2Nir + x 

V 0 sma>t if 2Nn + x £t < t a + 2Nit + ir 

-F p exp[-(f-/ 0 -2^-7r) 2 /r 2 ] if t 0 + 2Na + it <, t < 2Nit + 2ic 



F 0 sinfi>r 

V p exp[-(t-t 0 -2Nit)*/T 2 ] 



10 



In all of which equations El, E2 and E3, the parameters are to be interpreted as 
described herein. 



